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Limit-cycle  oscillators  are  used  to  model  a  broad  range  of  periodic  nonlinear  phenomena.  Using 
the  optically-injected  semiconductor  laser  as  a  paradigmatic  example,  we  demonstrate  that  at  spe¬ 
cific  operating  points,  the  period-one  oscillation  frequency  is  simultaneously  insensitive  to  multiple 
perturbation  sources.  In  our  system  these  include  the  temperature  fluctuations  experienced  by  the 
master  and  slave  lasers  as  well  as  fluctuations  in  the  bias  current  applied  to  the  slave  laser.  Tuning 
of  the  oscillation  frequency  then  depends  only  on  the  injected  optical  field  amplitude.  Experimen¬ 
tal  measurements  are  in  detailed  quantitative  agreement  with  numerical  modeling.  These  special 
operating  points  should  prove  valuable  for  developing  ultrastable  nonlinear  oscillators,  such  as  a 
narrow-linewidth,  frequency-tunable  photonic  microwave  oscillator. 


Nonlinear  dynamics  in  oscillators  has  been  invoked  to 
explain  a  wide  variety  of  physical  phenomena  over  disci¬ 
plines  ranging  from  neuroscience  to  geoscience  [1].  Sys¬ 
tems  and  devices  exhibiting  self-sustained  oscillations, 
described  mathematically  as  a  limit  cycle,  are  funda¬ 
mental  components  in  complex  systems  such  as  biolog¬ 
ical  oscillators  and  technological  applications  such  as 
time/frequency  references,  and  there  currently  is  consid¬ 
erable  interest  in  photonic  implementations  [2].  In  all  of 
these  instances,  there  is  the  need  to  understand  the  influ¬ 
ence  of  a  noisy/perturbing  environment  on  the  properties 
of  the  limit  cycle.  Perturbations  in  oscillatory  systems 
can  be  translated  in  frequency  through  the  nonlinear  cou¬ 
pling  among  elements.  For  example,  low-frequency  vibra¬ 
tion  or  temperature  perturbations  can  negatively  impact 
the  stability  of  a  high-frequency  oscillator. 

Here  we  demonstrate  that  nonlinear  dynamics,  which 
usually  degrades  system  performance,  can  counterintu¬ 
itively  be  used  to  suppress  the  deleterious  effects  such 
perturbations  have  on  performance  through  a  proper 
choice  of  the  operating  point  of  the  system.  Using  an 
optically-injected  semiconductor  laser,  we  confirm  that 
at  certain  operating  points  the  fundamental  resonance 
frequency  of  the  system,  the  so-called  period-one  (PI) 
frequency,  possess  greatly  reduced  sensitivity  to  current 
fluctuations  in  the  slave  laser  and/or  to  perturbations  in 
the  operating  temperatures  of  either  laser.  These  vari¬ 
ations  would  normally  lead  to  optical  frequency  fluctu¬ 
ations  of  the  free-running  lasers.  Effectively,  at  these 
points  the  tuning  of  the  PI  frequency  is  controlled  by  the 
amplitude  of  the  injected  optical  signal  alone.  This  is  in 
stark  contrast  to  the  beat  signal  generated  at  a  photo¬ 
diode  simultaneously  subject  to  two  free-running  lasers, 
where  the  frequency  of  the  current  oscillations  due  to  in¬ 
terference  is  controlled  simply  by  their  frequency  detun¬ 
ing.  This  system  provides  a  concrete  physical  platform  to 
investigate  the  role  of  nonlinear  dynamics  in  controlling 
sensitivity  to  external  perturbations,  and  should  there¬ 
fore  provide  an  additional  avenue  to  be  exploited  in  the 
creation  of  low  phase-noise  sources  e.g.  tunable  photonic 


oscillators. 

Model  systems  for  nonlinear  dynamics,  where  detailed 
quantitative  comparisons  between  model  and  the  real 
physical  system  can  be  made,  are  particularly  useful  to 
understand  the  complexities  introduced  by  nonlinearities 
in  the  actual  system.  Optically-injected  semiconductor 
lasers  have  been  shown  to  exhibit  a  wide  range  of  generic 
dynamic  characteristics  due  to  the  nonlinearities  inherent 
in  the  coupled  optical-field/free-carrier  system  [3].  These 
characteristics  have  been  quantitatively  recovered  using 
a  three-dimensional,  lumped-element,  detuned-oscillator 
model  of  the  coupling  between  the  circulating  optical  field 
and  the  carrier  density  in  the  gain  medium  [3-5]: 


X  = 

zx  +  ( bz  —  u)y  +  £ 

(1) 

y  = 

zy  —  (bz  —  uj)x 

(2) 

z  = 

k  —  Az  —  B(  1  +  2  z)(x2  +  y2  —  1) 

(3) 

These  equations  describe  a  system  of  three  coupled  ele¬ 
ments,  two  (x,y)  with  conjugate  coupling,  and  a  third, 
z,  exciting  the  other  two  while  being  damped  by  them, 
along  with  various  external  control  and  coupling/biasing 
parameters.  For  the  specific  case  of  the  optically-injected 
semiconductor  laser,  x  and  y  are  proportional  to  the 
quadrature  components  of  the  circulating  optical  field, 
and  z  is  proportional  to  the  carrier  density  of  the  gain 
medium.  The  parameters  £  and  w  describe  the  normal¬ 
ized  amplitude  and  detuning  (relative  to  the  free-running 
slave  laser)  of  the  injected  optical  field,  respectively,  while 
k  is  an  externally  applied  injection  current  beyond  the 
steady-state  bias.  The  parameter  b  is  intrinsic  to  the  in¬ 
jected  laser  and  is  the  ratio  of  the  real  and  imaginary 
parts  of  the  complex  refractive  index  variation  with  car¬ 
rier  density,  while  A  and  B  describe  coupling  terms  that 
depend  on  the  steady-state  bias  current.  All  parameters 
except  b  can  be  widely  varied  for  a  single  laser  experimen¬ 
tally  simply  by  changing  its  applied  current  and  operat¬ 
ing  temperature.  Nevertheless,  due  to  fabrication  issues, 
b  is  found  to  vary  between  lasers.  Therefore,  this  system 
is  capable  of  simulating  a  wide  range  of  external  coupling 
and  internal  cross-coupling.  The  key  control  parameters 
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FIG.  1:  Schematic  of  a  semiconductor  laser  subject  to  ex¬ 
ternal  optical  injection  (doted  green).  The  master  laser  is 
isolated  from  feedback.  In  our  experiment,  the  output  from 
the  injected  laser  (solid  red)  is  detected  by  a  fast  photodi¬ 
ode  and  monitored  using  a  microwave  spectrum  analyzer.  A 
small-signal  modulation  current  is  added  to  either  the  master 
or  slave  laser.  A  tunable  laser  acting  as  a  local  oscillator  (LO) 
can  be  added  to  generate  a  high-resolution  optical  spectrum 
using  a  heterodyne  technique  (dashed  blue). 


for  these  dynamics  are  the  amplitude  and  detuning  be¬ 
tween  the  injected  signal  from  the  master  laser  relative  to 
the  free-running  slave  laser  characteristics.  Over  a  wide 
range  of  these  two  control  parameters,  the  output  of  the 
injected  slave  laser  is  a  periodic  oscillation  at  the  PI  fre¬ 
quency  of  the  system  [3,  4]. 

Figure  1  is  a  schematic  of  the  configuration  of  an 
externally-injected  semiconductor  laser.  The  output  of  a 
steady-state  single-mode  master  laser,  at  an  optical  fre¬ 
quency  Vi ,  passes  through  an  isolator  and  is  injected  into 
the  oscillating  mode  of  a  second  semiconductor  laser  (the 
so-called  slave  laser),  also  steady-state  and  single  mode, 
with  a  free-running  frequency  v$.  A  weak  modulation 
current  can  be  added  to  either  the  master  or  slave  laser 
for  diagnostic  purposes.  The  slave  laser  under  investiga¬ 
tion  is  a  fiber-pigtailed  DFB  laser  operating  at  a  wave¬ 
length  of  1557  nm  that  has  been  extensively  investigated 
previously  [3] .  The  different  components  are  all  fiber  cou¬ 
pled  with  SMF-28  optical  fiber.  An  amplified  fast  photo¬ 
diode  detects  the  output  and  the  power  spectrum  of  the 
photodiode  signal  is  monitored  using  a  microwave  spec¬ 
trum  analyzer  (MSA) .  A  third  laser  is  sometimes  used  as 
a  tunable  local  oscillator  to  generate  a  heterodyne  signal 
that  effectively  converts  the  MSA  into  an  optical  spec¬ 
trum  analyzer  [3]. 

When  displaying  a  PI  dynamic,  the  optical  spectrum 
consists  of  dominant  peaks  at  the  injection  frequency 
and  at  one  resonance  frequency  (/o)  below,  which  is 
the  injection-perturbed  oscillation  frequency  of  the  slave 
laser.  In  addition,  other  weaker  peaks  at  periodic  offsets 
of  the  PI  frequency  are  typically  present.  The  addition 
of  a  weak  modulation  current  at  a  frequency  fm  <C  /o 
causes  sidebands  to  appear  about  these  main  optical 
peaks.  Figure  2  shows  typical  power  spectra.  The  power 
spectrum  displays  a  single  dominant  peak  at  the  PI  os¬ 
cillation/pulsation  frequency,  along  with  sidebands  and 
a  peak  at  the  modulation  frequency  if  a  modulation  cur¬ 
rent  is  present.  With  no  modulation  current,  only  the 
dominant  peak  at  the  PI  frequency  appears  in  the  power 
spectrum,  and  the  optical  spectrum  displays  only  the  dis- 


FIG.  2:  Power  spectra  of  the  laser  output  in  the  PI  oscilla¬ 
tion  regime  with  the  addition  of  a  weak  modulation  current 
to  the  slave  laser.  The  sidebands  are  minimized  when  the 
master  laser  is  detuned  to  generate  a  local  minimum  of  the 
PI  frequency  (solid  red,  detuning  =  -2.1  GHz),  but  are  much 
stronger  for  a  shifted  detuning  frequency  (dashed  blue,  detun¬ 
ing  =  -1.2  GHz).  A  normalized  injection  amplitude  of  0.06  is 
used  in  both  cases  and  the  peak  at  the  modulation  frequency 
(500  MHz)  remains  essentially  unchanged. 


Crete  set  of  peaks  offset  by  this  frequency. 

For  these  experiments,  the  operating  points  of  the  mas¬ 
ter  and  slave  lasers  are  temperature  tuned  so  the  PI  fre¬ 
quency  is  locally  insensitive  to  the  master-slave  detun¬ 
ing;  the  PI  frequency  is  at  a  local  minimum  with  respect 
to  this  control  parameter  and,  therefore,  insensitive  to 
temperature  fluctuations.  A  weak  sinusoidal  modulation 
current  with  frequency  fm  <C  /o  is  then  added  to  either 
the  master  or  slave  laser.  Typically,  the  PI  frequency 
peak  develops  sidebands  offset  by  integer  multiples  of 
fm  regardless  of  which  laser  the  modulation  current  is 
added  to.  However,  at  approximately  the  detuning  that 
produces  the  minimum  PI  frequency  we  observe  a  mini¬ 
mum  in  the  amplitude  of  the  modulation  sidebands  in  the 
power  spectrum  when  the  modulation  current  is  added 
to  the  slave  laser,  as  shown  in  Fig.  2.  By  comparison,  the 
peak  at  the  modulation  frequency  shows  no  such  mini¬ 
mum,  verifying  that  the  linear  response  of  the  laser  is 
unchanged.  If  the  modulation  current  is  added  to  the 
master  laser,  there  is  a  less  pronounced  modulation  min¬ 
imum,  but  it  occurs  at  a  detuning  well  offset  from  the 
operating  point  where  the  PI  frequency  is  minimized. 

Figure  3  shows  experimental  measurements  of  the  am¬ 
plitude  of  the  sidebands  and  the  PI  frequency  as  the  de¬ 
tuning  frequency  of  the  master  laser  is  stepped  through 
the  PI  frequency  minimum.  Here  the  amplitude  of  the 
PI  frequency  peak  varies  strongly  as  the  detuning  moves 
away  from  the  Hopf  bifurcation.  Therefore,  the  ampli¬ 
tudes  of  the  sidebands  are  shown  normalized  to  the  am¬ 
plitude  of  the  central  peak  at  each  detuning.  The  error 
for  the  experimentally  measured  detuning  is  ±300  MHz 
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Detuning  Frequency  (GHz) 


FIG.  3:  Experimentally  measured  (symbols)  and  calculated 
(solid  curves)  values  of  the  modulation  sideband  amplitude 
and  the  PI  frequency  using  a  normalized  injection  ampli¬ 
tude  of  0.06.  This  figure  shows  the  positive  frequency  de¬ 
tuning  sideband  (red  circles),  the  negative  detuning  sideband 
(green  diamonds),  and  the  corresponding  PI  frequency  (blue 
squares) . 


while  the  depths  of  the  minima  were  limited  by  the  back¬ 
ground  noise  levels.  The  dips  in  the  amplitudes  of  the 
sidebands  are  quite  narrow  with  respect  to  the  relaxation 
rates  of  the  free-running  laser  which  are  >  109s-1  at  this 
operating  point  [3,  5].  Further,  additional  measurements 
made  on  signals  of  broader  bandwidth  (>  10  MHz)  up  to 
modulation  frequencies  of  100  MHz  showed  an  equivalent 
dip  in  the  amplitude.  Therefore,  there  is  reduced  sen¬ 
sitivity  to  all  low-frequency  fluctuations  in  the  injection 
current  of  the  slave  laser  simultaneously  with  the  reduced 
sensitivity  to  fluctuations  in  the  master-slave  detuning. 

Figure  4  shows  a  mapping  of  key  operating  points  of 
the  laser,  comparing  experimental  measurements  with 
model  calculations  (described  later),  as  functions  of  the 
master  laser  detuning  from  the  free-running  slave  optical 
frequency  and  the  amplitude  of  the  injected  signal.  The 
mapping  is  made  by  monitoring  the  power  and  optical 
spectra  of  the  optically-injected  laser.  Saddle-node  and 
Hopf  bifurcations  bound  the  region  of  stable  locking,  with 
the  Hopf  bifurcation  leading  to  PI  limit-cycle  oscillation. 
More  complex  dynamics  occur  in  the  region  between  the 
points  labeled  as  either  the  saddle-node  bifurcation  of 
limit  cycles  or  alternate  routes  to  chaos  (through  period 
doubling  or  other  routes).  One  cannot  identify  the  sad¬ 
dle  node  bifurcation  of  limit  cycles  through  the  spectra, 
but  previous  work  using  bifurcation  analysis  has  done  so 
[4] .  Along  the  other  boundary,  period-doubling  routes  to 
chaos  have  been  identified  [3]. 

To  produce  the  quantitative  comparisons  between 
data  and  theory  we  modeled  the  optically-injected  laser 
using  a  set  of  coupled  equations,  of  the  form  of  Eqs.  1-3, 
for  the  circulating  field  amplitude  of  the  oscillating  laser 
mode  and  the  carrier  density  of  the  gain  medium  [3,  6]. 


FIG.  4:  Mapping  of  key  transitions  of  the  injected  laser  as  a 
function  of  the  detuning  frequency  and  amplitude  of  the  mas¬ 
ter.  Experimental  Data  -  the  open  squares  are  the  Saddle- 
Node  bifurcation  points  that,  along  with  the  green  triangles 
representing  Hopf  Bifurcation  points,  bound  the  region  of  sta¬ 
ble  locking,  the  red  diamonds  mark  operating  points  with  a  lo¬ 
cal  PI  minimum  frequency.  The  pink  circles  denote  a  saddle- 
node  bifurcation  of  limit  cycles,  either  from  PI  to  more  com¬ 
plex  dynamics  or  between  two  PI  frequencies,  and  the  purple 
diamonds,  alternate  routes  from  PI  to  more  complex  dynam¬ 
ics.  Calculated  Data  -  the  light  blue  regions  mark  PI,  and 
the  dark  blue  regions  P2,  periodic  dynamics  with  lines  of  con¬ 
stant  frequency  are  labeled  in  GHz.  The  white  region  denotes 
stable  operation,  while  the  black  represents  complex  dynam¬ 
ics.  The  black  squares  denote  the  calculated  positions  of  the 
local  PI  minima  to  be  compared  with  the  red  diamonds. 


This  model  modifies  Eqs.  1-3  by  adding  the  effects  of  gain 
saturation  to  both  the  real  (refractive  index)  and  imagi¬ 
nary  (gain)  parts  of  the  nonlinear  susceptibility  [3] .  The 
numerical  simulations  used  a  Runge-Kutta  integration 
procedure  over  durations  >  1/rs  to  generate  time  series 
that  were  subsequently  Fourier  transformed  to  produce 
spectra  [6]. 

The  results  of  the  numerical  calculations  are  shown  in 
Figs.  3  and  4  using  parameters  that  have  been  determined 
experimentally  for  the  laser  under  investigation  [5].  The 
model  fully  reproduces  the  bifurcation  boundaries  of  the 
PI  region,  the  local  minimum  in  the  PI  frequency,  and 
the  simultaneous  dip  in  the  amplitude  of  the  sidebands  as 
the  detuning  between  master  and  slave  lasers  is  changed. 
The  positions  of  the  PI  frequency  minimum  and  the  am¬ 
plitude  dip  are  sensitive  to  the  effects  of  gain  satura¬ 
tion;  to  obtain  good  agreement  between  experimentally 
measured  data  and  the  model  saturation  effects  were  in¬ 
cluded.  More  details  on  the  comparison  of  the  model 
and  experimental  data,  particularly  optical  spectra,  will 
be  given  in  a  later  publication.  Here,  we  note  that  the 
basic  observation  of  simultaneous  insensitivity  to  pertur¬ 
bations  in  the  detuning  between  master  and  slave  and  to 
perturbations  in  the  slave  laser  current  remained  when 


3 

Approved  for  public  release;  distribution  unlimited. 


4 


gain  saturation  effects  were  removed,  confirming  that  this 
term  in  the  model  is  not  necessary  to  generate  the  under¬ 
lying  interaction,  only  to  achieve  quantitative  agreement 
between  our  numerical  model  and  experiment. 

It  is  important  to  remember  that  the  sidebands  in  the 
power  spectrum  are  the  sum  of  terms  due  to  the  mix¬ 
ing  of  the  sidebands  in  the  optical  spectrum  with  ad¬ 
jacent  central  peaks.  By  comparing  the  power  spectra 
with  optical  spectra  generated  simultaneously  using  the 
hetrodyne  technique,  we  verified  that  the  sharp  minima 
in  the  sidebands  (see  Fig.  3)  result  primarily  from  the 
interference  of  the  multiple  wave-mixing  contributions, 
not  the  simultaneous  reduction  in  the  amplitude  of  the 
sidebands  of  each  optical  peak.  The  component  at  the 
modulation  frequency,  which  results  from  the  sum  of  the 
mixing  of  each  optical  sideband  with  its  respective  cen¬ 
tral  peak,  remains  essentially  constant  as  the  injection 
detuning  frequency  is  changed. 

While  the  effects  of  current  noise  on  the  PI  frequency 
are  suppressed,  spontaneous  emission  noise  still  plays  a 
role.  One  can  see  this  in  the  model  by  considering  spon¬ 
taneous  emission  as  a  fluctuating  source  term  with  con¬ 
tributions  to  both  the  optical  field  and  carrier  density 
equations  [7].  The  carrier-equation  contribution  is  elimi¬ 
nated  through  the  interference  of  wave-mixing  terms  de¬ 
scribed  above,  as  if  it  were  an  injection  current  modula¬ 
tion  source,  while  the  field-equation  contribution,  which 
acts  like  a  stochastic  external  optical  injection  term,  still 
contributes.  However,  only  the  amplitude  fluctuations, 
not  the  phase/frequency  fluctuations,  contribute,  lead¬ 
ing  to  the  Lorentzian  lineshape  of  the  spectral  features  in 
Fig.  2.  Previously,  it  was  shown  how  the  coupling  of  opti¬ 
cal  fields  and  free  carriers  in  a  semiconductor  laser  could 
lead  to  reduced  power  fluctuations  [7].  In  the  optically- 
injected  laser,  the  nonlinear  dynamics  can  completely  flip 
the  sensitivity,  so  that  it  is  the  frequency  fluctuations 
that  are  minimized. 

The  curve  of  measured  and  modeled  operating  points 
in  Fig.  4  where  the  PI  frequency  is  minimized  is  initi¬ 
ated  along  a  saddle-node  bifurcation  of  limit  cycles.  In 
fact,  at  the  end  of  the  bifurcation  the  transition  is  be¬ 
tween  two  limit  cycles  of  different  frequency.  The  laser 
switches  from  a  limit  cycle  where  the  injected-shifted 
laser  oscillation  frequency  is  pushed  away  from  the  injec¬ 
tion  frequency  at  larger  detuning  to  a  limit  cycle  where 
the  oscillation  frequency  is  pulled  towards  the  injection 
frequency  for  smaller  detuning.  This  change  in  relative 
behavior  of  the  system  resonance  relative  to  the  injection 
frequency  is  necessary,  though  at  higher  injection  levels 
we  observe  a  smooth  rather  than  abrupt  transition.  Fur¬ 
ther,  we  observe  that  the  minimum  P 1  frequency  merges 
with  the  bifurcation  line  as  the  difference  between  the 
two  frequencies  goes  to  zero.  Frequency  pushing  behav¬ 
ior  arises  from  the  fact  that  the  semiconductor  laser  is  a 
detuned  oscillator,  while  the  frequency  pulling  behavior 
is  the  kind  expected  from  an  Adler-type  analysis  of  an  in¬ 
jected  oscillator  [8,  9].  At  low  injection  levels,  the  saddle- 


node  of  limit  cycle  bifurcation  starts  at  a  PI  frequency 
and  master  laser  frequency  offset  approximately  equal  to 
the  characteristic  relaxation  resonance  frequency  of  the 
free-running  laser,  separating  the  frequency-pushing  and 
pulling  regimes.  At  higher  injection  levels,  where  the  sen¬ 
sitivity  to  perturbations  is  reduced,  the  local  minimum 
of  the  PI  frequency  is  also  a  signature  of  a  localized  re¬ 
gion  of  relative  frequency  pulling,  or  reduced  absolute 
frequency  pushing,  of  the  injection-perturbed  laser  oscil¬ 
lation  frequency. 

At  the  high-injection  level  end  of  the  operating  points 
where  the  PI  frequency  is  minimized,  the  Hopf  bifur¬ 
cation  is  crossed  at  approximately  the  position  where 
it  deviates  most  strongly  from  the  line  ui  =  £  that 
the  bifurcation  asymptotically  approaches  with  increas¬ 
ing  injection  [10].  Therefore,  we  observe  that  the  re¬ 
gion  of  reduced  sensitivity  operation  is  bounded  within 
a  range  where  the  attractor  characteristics  are  different 
(frequency-pulling)  from  the  PI  characteristics  at  larger 
detuning  and  away  from  the  regions  of  complex  dynam¬ 
ics  and  the  saddle-node  bifurcation  of  limit  cycles  (fre¬ 
quency  pushing).  Analysis  of  the  PI  dynamics  in  terms 
of  a  solution  with  two  dominant  optical  frequency  com¬ 
ponents  [11]  cannot  yield  the  localized  minima  of  the  PI 
frequency  with  changes  in  the  detuning.  Similarly,  a  lin¬ 
earized  analysis  of  the  laser  under  stable  injection  lock¬ 
ing  [10]  cannot  yield  such  minima  in  this  region,  though 
they  appear  under  the  full  nonlinear  analysis  [12].  There¬ 
fore,  the  PI  region  and  even  the  stable  locking  region 
can  have  attractors  that  yield  dynamical  features  un¬ 
expected  from  perturbation  solutions.  Nonetheless,  the 
full  nonlinear  model,  Eqs.  1-3,  recovers  these  features 
qualitatively/semi-quantitatively,  with  a  proper  account¬ 
ing  of  the  gain/refractive  index  saturation  being  nec¬ 
essary  to  quantitatively  recover  them  in  the  optically- 
injected  semiconductor  laser. 

In  summary,  we  have  presented  experimental  mea¬ 
surements  and  numerical  calculations  based  on  the  cou¬ 
pled  optical- field/carrier-density  model  of  the  optically- 
injected  semiconductor  laser  that  show  the  existence  of 
specific  operating  points  with  reduced  sensitivity  to  sys¬ 
tematic  fluctuations.  These  techniques  complement  and 
expand  upon  the  recent  use  of  specifically  engineered 
nonlinear  oscillators  that  used  nonlinear  dynamics  to 
suppress  oscillator  phase  noise  [13,  14].  Further,  the  ex¬ 
cellent  agreement  between  model  and  experiment  in  this 
system  makes  it  an  ideal  test  configuration  for  investi¬ 
gating  novel  responses  of  a  nonlinear  system  to  exter¬ 
nal  stimulus.  More  generally,  the  work  highlights  the 
nontrivial  changes  in  the  response  of  a  nonlinear  system 
to  perturbations,  and  the  fact  that  at  specific  operating 
points  nonlinearly  shifted  perturbations  are  suppressed, 
while  otherwise  appearing  to  remain  similar  in  character¬ 
istics  to  nearby  operating  points.  This  has  technological 
relevance  to  the  frequency  reference  application  cited  ear¬ 
lier,  but  we  believe  that  such  localized  operating  points 
may  also  be  of  importance  to  a  wide  range  of  nonlinear 
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